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Abstract 



We apply a Langevin model to investigate the nonlinear ac responses 

of electro-magnetorheological (ERMR) fluids under the application of two 

crossed dc magnetic (z axis) and electric (x axis) fields and a probing ac 

sinusoidal magnetic field. We focus on the influence of the magnetic fields 

which can yield nonlinear behaviors inside the system due to the particles with 

a permanent magnetic dipole moment. Based on a perturbation approach, we 

extract the harmonics of the magnetic field and orientational magnetization 

analytically. To this end, we find that the harmonics are sensitive to the 

degree of anisotropy of the structure as well as the field frequency. Thus, it 

is possible to real-time monitor the structure transformation of ERMR fluids 

by detecting the nonlinear ac responses. 
PACS: 83.80.Gv, 72.20.Ht, 41.20.-q, 82.70.Dd 
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I. INTRODUCTION 



Electro-magnetorheo logical (ERMR) fluids [1,2] behave as both electrorheological (ER) 
fluids [3-6] and magnetorheological (MR) fluids [7-9], and are, in general, particle suspen- 
sions where the particles have large dielectric constants and permanent magnetic moments. 
In fact, there exist a variety of particles which can be polarized by both electric field and 
magnetic field [10]. For instance, there is one candidate for this system, namely, a suspension 
of titanium-coated iron particles in a silicon oil [11]. As the external electric field or mag- 
netic field exceeds a threshold, ER or MR fluids turn into a semi-solid, the ground state of 
which is a body-centered-tetragonal (bet) lattice. For an ERMR solid, one [1] proposed that 
a structure transformation of the ERMR solid from the bet lattice to some other lattices can 
appear when a magnetic field is simultaneously applied perpendicular to the electric field. 
To one's interest, such a structure transformation from bet to face-centered-cubic lattice was 
experimentally observed as the ratio between the magnetic and the electric fields exceeded 
a minimum value [2] . Recently, one of the present authors showed that an alternative struc- 
ture transformation from the bet structure to the fee can appear under the application of 
electric fields only [12]. 

Electric or magnetic fields of high strength applied to the ERMR solid produce a nonlin- 
earity in the dependence of the polarization or magnetization on the field strength. Conse- 
quently, in the presence of an ac electric or magnetic field, the electrical or magnetic response 
will in general consist of ac fields at frequency of the higher-order harmonics [13-16]. A con- 
venient method of probing the nonlinear characteristics of the composite is to measure the 
harmonics of the nonlinear polarization (or magnetization) under the application of a sinu- 
soidal ac field [17]. In this case, the strength of the nonlinear polarization or magnetization 
should be reflected in the magnitude of the harmonics. For extracting such harmonics, the 
perturbation approach [18,19] and self-consistent method [19,20] can be used. To the best 
of our knowledge, no work has been done on the nonlinear ac responses of ERMR fluids. 

For the present system under consideration, the nonlinearity can be caused to appear 



2 



by two effects, namely, normal saturation and anomalous saturation. In detail, the normal 
saturation arises from the higher terms of the Langevin function at large field intensities [21]. 
In contrast, the anomalous saturation results from the equilibrium between entities with 
higher and lower dipole moments which is shifted under the influence of the field [21]. 
Our formalism will hold for the coupling between the normal saturation and anomalous 
saturation. 

In the present paper, to investigate the structural effect on the nonlinear ac responses 
of ERMR solids, we shall apply a Langevin model to derive the orientational magnetization 
as well as the magnetic field inside the ERMR solid. In this connection, the perturbation 
approach will be used to to extract the harmonics of the magnetic field and orientational 
magnetization. 

This paper is organized as follows. In Sec. II, we apply the Langevin model to derive 
the orientational magnetization as well as the magnetic field inside the ERMR slid, and 
calculate the effective permeability of the ERMR solid by using a generalized Clausius- 
Mossotti equation. Also, based on the perturbation approach, we extract the harmonics of 
the magnetic field and orientational magnetization analytically. In Sec. Ill, we numerically 
investigate these harmonics as a function of the degree of anisotropy of the structure, as well 
as the frequency of the ac magnetic field. This is followed by a discussion and conclusion in 
Sec. IV. 

II. FORMALISM 

A. Nonlinear characteristics 

Under the application of a strong magnetic field along z axis, a nonlinear characteristic 
can appear in the ERMR fluid, due to the particles with a permanent magnetic dipole 
moment p . Accordingly, the dependence of the magnetic induction B on the field H will 
be nonlinear [21] 
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B = /x e H + 47rx#o H o, 



(1) 



where x an d A*e stand for the nonlinear susceptibility and effective permeability for the 
longitudinal field case, respectively. In this case, the effective permeability /i e is determined 
by the generalized Clausius-Mossotti equation [12] 



where \i2 represents the permeability of the host fluid, ~N\ the number density of the par- 
ticles, ks the Boltzmann constant, T the absolute temperature, / the frequency of the 
applied magnetic field, T\ the relaxation time of the particles, and a.\ the magnetizability 
of the particles. In Eq. (2), the longitudinal demagnetizing factor g L deserves a thorough 
consideration. For an isotropic array of magnetic dipoles, the demagnetizing factor will be 
diagonal with the diagonal element gi = 1/3. However, in an anisotropic array like ERMR 
solids, the demagnetizing factor can still be diagonal, but it deviates from 1/3. In fact, the 
degree of anisotropy of the system is just measured by how g L is deviated from 1/3. It is 
worth noting that gL < 1/3 in the present longitudinal field case. Furthermore, there is a 
sum rule for the factors, g L + 2g T = 1 [19,22], where gx denotes the transverse demagnetizing 
factor. Such factors were measured by means of computer simulations [23,24]. Thus, to in- 
vestigate the anisotropic structural information of the array, we have to modify the Clausius- 
Mossotti equation accordingly by including the demagnetizing factor. The substitution of 
9l(— 9t) = 1/3 into Eq. (2) yields the usual (isotropic) Clausius-Mossotti equation, which 
does not include the particle-particle interaction. In fact, when we studied the field-induced 
structure transformation in ER solids, we developed the generalized Clausius-Mossotti equa- 
tion [12] by introducing a local-field factor f3' which reflects the particle-particle interaction 
between the particles in a lattice. In detail, the generalized Clausius-Mossotti approach 
[Eq. (2)] is a self-consistent determination of the local field due to a lattice of dipole mo- 
ments. That is, Eq. (2) should be expected to include the particle-particle interaction, and 
the degree of the particle-particle interaction depends on how much g^ deviates from 1/3 
(note that gi = P'/S)- Indeed our numerical results will also show that as gL deviates from 
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1/3 more (namely, more and more particle chains are formed, and the particle-particle in- 
teraction becomes more and more strong, too), the obtained harmonics become more large 
accordingly, see Figs. 1-2. These harmonics just reflect the magnitude of the nonlinear 
susceptibility, as expected. 

In Eq. (2), the term Pq/(3/c b T) results from the average contribution of the permanent 
magnetic dipole moment to the average value of the work required to bring a particle into 
the field H . More precisely, the mean value of the component of the permanent dipole 
moment in the direction of the field is given by 



Regarding this linear relation in use, we should make some remarks. In the present work, 
we shall adopt the perturbation approach [25], which is suitable for weak nonlinearity. In 
this perturbation approach, it is well established that the effective third-order nonlinear 
susceptibility can be calculated from the linear field [26], while the effective higher-order 
nonlinearity must depend on the nonlinear field [25]. Alternatively, we could adopt the self- 
consistent method [19], but the perturbation approach appears to be more convenient for 
analytic expressions [19]. Thus, for focusing on (weak) third-order nonlinearity, it suffices 
to use the Clausius-Mossotti equation [Eq. (2)] by taking into account of the linear relation 
[Eq. (4)] only. Due to the same reason, the contribution from the nonlinear field will be 
omitted throughout the paper. 

In view of Eq. (1), the desired field-dependent incremental permeability /in is given 




(3) 



with 7 = p H /(k B T). That is, we set the Langevin function 



L( 7 ) = 7 /3. 



(4) 



by [21] 



(5) 



Then, the nonlinear magnetic effect is characterized by A/j/Hq [21] 
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= 0-2 = 127r X- (6) 

In the present work, the magnetization can be split up to two parts, i.e., the induced 
magnetization due to the magnetizability of the particles and the orientation magnetization 
M or due to the alignment of permanent magnetic dipole moments inside the particles. Let 
us now consider the orientation magnetization M or of a sphere of volume V, containing 
rii particles with a permanent magnetic dipole moment [see Eq. (8)] embedded in a 
continuum with permeability ji^ [see Eq. (9)]. The sphere is surrounded by an infinite 
medium with the same macroscopic properties as the sphere. 

The average component in the direction of the field of the magnetic dipole moment due 
the the dipoles in the sphere (M d • e) = VM or is given by: 

{ m d . e) = V M„ = > AX !A ^ e e : ikBT , (7) 

where e denotes the unit vector in the direction of the external field, and X stands for 
the set of position and orientation variables of all particles. Here u is the energy related 
to the dipoles in the sphere, and it consists of three parts: the energy of the dipoles in 
the external field Ud e , the magnetostatic interaction energy of the dipoles u mi , the non- 
magnetostatic interaction energy between the dipoles u nm i which is responsible for the short- 
range correlation between orientations and positions of the dipoles. Regarding the particle- 
particle interaction energy, Eq. (2) should be expected to contain both u mi and u nmi (at 
least to some extent) as g L is not equal to 1/3. Nevertheless, for the numerical calculations 
in the next section, Eqs. (39) and (40) will be used which mean that u nmi is omitted due to 
the predominant external magnetic field. In Eq. (7), is given by 

ni 

M d = J2(Pdh (8) 

i=i 

with p d = po(/ioo + 2/i2)/3/i2, where /ioo represents the permeability at frequencies at which 
the permanent dipoles cannot follow the changes of the field but where the atomic and 
the electronic magnetization are still the same as in the static field. Therefore, fj,^ is the 
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permeability characteristic for the induced magnetization. In practice, //qo can be expressed 
in the expression containing an intrinsic dispersion, 

= "-<«» + TWJZ' (9) 

where /ioo(0) is the high-frequency limit permeability, and A/i stands for the magnetic dis- 
persion strength with a characteristic frequency f c . 

The external field in this model is equal to the field within a spherical cavity filled with 
a continuum of permeability fi^, while the cavity is situated in a medium with permeability 
fi e . This field is called Frohlich field Hp, as given by the Frohlich model [21,27] 

In this equation, higher-order terms have been omitted. Regarding the derivation of the 
nonlinear Frohlich field [Eq. (10)], here we would like to add some comments. For the 
spherical cavity under consideration, the linear cavity field (i.e., linear Frohlich field) H F n 
can be easily obtained by solving the usual electrostatic equation. That is, it should be 
the first term of the right-hand side of Eq. (10). Next, to obtain the present nonlinear 
Frolich field H F , following Ref. [21] we took one step forward to develop H F n in a Taylor 
series around fi e and neglecting terms in Hq and higher powers of H . It is worth noting 
that the higher powers of H can also contribute to the resulting harmonics. However, this 
contribution is small enough to be neglected. 
Now we use 

°" H/ <>• (11) 



dH F 
Then, we obtain 

(M d - e)\ HF=0 = 0, (12) 
8 <M d ■ e) = 7^[((M d ■ e) 2 ) - (M d • e) 2 ], (13) 



d HF s k B T [ 

^(M,.e)U F=0 = -l-(M 2 ) , (14) 



<9 3 1 

■{M d ■ e) = — — -[{(M d • e) 4 > - 3((M d • e) 3 ) 



dH F x a 1 (k B Tf [ 

(M d ■ e) + 6((M d • e) 2 )(M d • e) 2 - 3((M d • e) 2 ) 2 



6(M d • e) 4 + 6(M d ■ e) 2 ((M d • e) 2 ) - (M d ■ e)((M d ■ e) 3 )], (15) 
d 3 1 

dHf {Md ' e)lHp=o = mk~ B Tr 



-(M, • e)\ Hr=0 = —r ^J 3(M 4 ) - 5(M 2 ) 2 ]. (16) 



Note the subscript indicates the absence of the field. 

To express the macroscopic saturation behavior in terms of microscopic quantities, one 
should take into account the higher derivatives of the average moment such that 

d(M d ■ e) | l d 3 (M d .e) 3 

(M d ■ e) = — \ HF=0 H F + - ^ \ Hf=0 H f 

li e (M 2 ) 47TX/ioo (M 2 ) ^3 
-"0 + 77 i 7t~, — 77T n ( 



2/i e + /ioo k B T (2/i e + /i^) 2 k B T 

27/x 3 3(M 4 ) - 5(il 
(2/i e + z^) 3 90(A; B T) 3 



o 



27/x 3 3(M 4 ) -5(M 2 ) 2 



where higher-order terms than 3rd have been neglected. 
In addition, we have a general relation 



{M d ■ e) = VM or = ^J^VH + X VHl (18) 

47T 



In view of Eq. (6) and the terms in Hq and H , we obtain 
Afi 18tt /i 4 



H 2 5(fc B r) s (2^ + |ioo)W + ^) 
On the other hand, we may write 



[3(«>o/VO - 5«JW^/V)]. (19) 



(^./^y^^^,, (20) 



(-^i)o /a'oo + 2/i 2 \ 4 rWi 4p/ a a \i /o-n 

— r^= — ^ [t/PoL( cos %Z^ cos ^ s >]. (21) 



In the right-hand sides of both Eqs. (20) and (21), the (• • -)'s denote 

/f rl V n i-« . . . p -u/k B T 

Hence the nonlinear magnetic increment A/i/Hq [Eq. (19)] can be explicitly expressed. 



B. Nonlinear magnetization and high-order harmonics 



1. Longitudinal field 

For the longitudinal field case, there is H (t) = H dc + H &c (t) = H dc + if ac sin cut along z 
axis, with u = 2nf '. Here H dc denotes the dc field which induces the anisotropic structure 
in the ERMR solid, and i? ac (t) stands for a sinusoidal ac field. Then, the Frohlich field is a 



superposition of odd- and even- order harmonics such that 

H F = H { p c) + H u sin ujt + E 2uJ cos 2ut + H 3u> sin 3cut + • • • . (23) 

Accordingly, the orientational magnetization contains harmonics as 

M or = M^ c) + M w sin ut + M 2u] cos 2ut + M 3uJ sin 3ut -\ . (24) 

where the harmonics (H u , H 2lu , H 3u ,, M u , M 2uJ , and M 3uJ ) and the dc components (H^ and 
M^) can be expressed as 

= H dc h + *-Hl c H dc h + Hlh, (25) 

H u = H^h + \Hlh + ZH^Hih, (26) 
3 

H 2lu = —-H^Hfoh, (27) 

= --H^ C I 3 , (28) 

= H dc J, + \HlH dcX + Hlx, (29) 

M w = H^J, + ^Hix + 3H, c H 2 dcX , (30) 

M 2w = -^if a 2 c iJ dcX , (31) 

M 3w = -^ a 3 c x, (32) 

with 7i = 3// e /(2/i e + /loo), J 3 = 127rx/ioo/(2/i e + /ioo) 2 , and Ji = (l/47r)(/i e - /i^). In 



the above derivation, we used two identities, i.e. sin 2 cut = (1 — cos2cut)/2 and sin 3 cut = 
(3/4) sin cut - (1/4) sin3cut 
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2. Transverse field 



For the transverse field case, only a sinusoidal electric field H (t) = sin cot is applied 
along x axis. In this connection, the Frohlich field is a superposition of odd-order harmonics 

H F = sin cot + H 3u> sin 3ut + • • • . (33) 

Accordingly, the orientational magnetization contains harmonics as 

M or = M w sin ut + M 3w sin 3ujt H . (34) 

In Eqs. (33) and (34), the harmonics H u , H 3u) , M w , M 3ul can be expressed as 

H„ = H & J[ + ^/ 3 , (35) 

^ = -^aVs, (36) 
M w = tf ac J( + ^Hlx', (37) 
M 3w = -\Hlx'- (38) 

with = 3/4/(2/4 + /^), 7 3 = 127rx'/ioo/(2/i^ + / u 00 ) 2 , and J[ = (l/47r)(/i' e - /i^), where // e 
and x' denote the effective permeability and nonlinear susceptibility in the transverse field 
case. For the transverse field case, x! is induced to appear by the moderate ac magnetic field 
only. Thus, in view of Eqs. (7) and (8), we find x! <^ X (x i s caused to appear by both the 
moderate ac and high dc magnetic fields in the longitudinal field case, as mentioned above). 
As a result, the third-order harmonics of the Frohlich field and orientational magnetization 
are much smaller in the transverse field case than in the longitudinal field case. In this 
regard, for the following numerical calculations, we shall focus on the longitudinal field case 
only. 

III. NUMERICAL RESULTS 

Now we are in a position to do some numerical calculations. For the numerical calcula- 
tion, we use 
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H(cos%) = 1, 



(39) 




^(cos%^^cos6» rs ) = -(5ni -2). 



(40) 



j=l r=l s=l 



Eqs. (39) and (40) imply that the correlation between the dipole moments of the particles 
and in turn the anomalous saturation are neglected in the sense that the current anomalous 
saturation is much more weak. This is reasonable because, in the presence of a magnetic 
field, the permanent magnetic moments of the particles are easily directed along the field. 
Therefore, the equilibrium between the higher dipole moment (of particle chains) and the 
lower dipole moment (of particle chains) might not be able to predict significant nonlinearity, 
when compared to the normal saturation. 

Next, we take one step forward to obtain the nonlinear magnetic increment 



For numerical calculations, take the following parameters: H^ c = 100 Oe, if ac = 1 Oe, 
//oo(O) = 1.1, A/i = 8, f c = 4.8 x 10 3 Hz, // 2 = 1, Po = 10" 13 emu, n = 4.8 x 10~ 7 s, and 
«i = 10~ 10 cm~ 3 . In addition, the radius of particle is taken to be 5 /im, and the volume 
fraction of the particles is 0.2. 

Fig. 1 shows the fundamental, second- and third- order harmonics of the Frohlich field as 
a function of the field frequency, for various gi in the longitudinal field case. In this figure, a 
peak is observed always due to the existence of an intrinsic dispersion [Eq. (9)]. In particular, 
as the system changes from isotropic case (ql = 1/3) to anisotropic (ql ^ 1/3) because of 
the appearence of the particle chains, the harmonics of the field can be changed accordingly. 
In detail, stronger anisotropy (namely, decreasing the longitudinal demagnizing factor g£) 
leads to larger harmonics, especially in the low-frequency region. Similar effect can be shown 
in Fig. 2 where the harmonics of the orientational magnetization in the longitudinal field case 
are investigated for various gi as well. However, the fundamental harmonic of the orientation 
magnetization [see Fig. 2(a)] behaves in a different way from that of the Frohlich field [see 




Hi ~ 5(k B Tr (2^ + ^)2(2^ + ^) 




(41) 
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Fig. 1(a)]. In detail, as the frequency increases, this fundamental harmonic decreases first, 
then increases, and after reaching a peak they decreases again. 

For the transverse field case, it could be concluded that as the demagtizing factor ql 
decreases, both the fundamental and third-order harmonics of the Frohlich field and orien- 
tation magnetization are caused to decrease accordingly, which is just opposite with those 
obtained from the longitudinal field case (no figures shown here). The reason is that there 
is a sum rule between gi and gr, gi + 2gr = 1- 

In a word, for the longitudinal field case, besides the odd-order harmonics, the even-order 
harmonics are also induced to appear due to the coupling between the applied dc and ac 
magnetic fields along z axis [see Figs. 1 and 2, or Eqs. (23)~(32)] even though only the cubic 
nonlinearity [Eq. (1)] is considered due to the virtue of symmetry of the system. On the 
other hand, for the transverse field case (i.e., there is only a single ac magnetic field applied 
along x axis), in view of the cubic nonlinearity [Eq. (1)] of interest, only the odd-order 
harmonics are induced, as already predicted by Eqs. (33)~(38). 

In addition, even though there is no particle-particle interaction (i.e., u mi = u nmi = 0) 
as gL = 1/3, the nonlinear behavior due to the normal saturation could still be induced to 
occur because of the presence of external fields, i.e., u de ^ 0. This is the reason why the 
harmonics shown in Figs. 1-2 are nonzero at gL — 1/3. 

From Figs. 1 and 2, we find that the second-order harmonics of the Frohlich field and 
orientational magnetization are of three orders of magnitude larger than the corresponding 
third-order harmonics. Thus, to monitor the structure transformation of ERMR solids, it is 
more attractive to detect the second-order harmonics than the third-order. 

Finally, we display the temperature effect on the harmonics of the Frohlich field and 
orientational magnetization in Fig. 3 and Fig. 4, respectively. It is shown that decreasing 
the temperature T causes all the harmonics to increase because of the change in the Frohlich 
field. 
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IV. DISCUSSION AND CONCLUSION 



Here some comments are in order. Based on a Langevin model, we have investigated the 
nonlinear ac responses of ERMR solids which are subjected to a structural transformation. 
In the present work, we have focused on the effect of the magnetic fields on the nonlinear 
ac responses. In fact, for the ERMR solid, the influence of electric fields can be investigated 
as well. However, the present Langevin model is invalid for this purpose because there 
is usually no permanent electric dipole moment inside the suspended particles, except for 
ferroelectric particles. To one's interest, when the particles own a nonlinear characteristic 
inside them, the electric-field effect on the nonlinear ac responses can still be discussed, 
based on the perturbation approach and self-consistent method [19] accompanying with an 
effective medium theory like the Maxwell- Garnett approximation. 

In this paper, the Frohlich field is actually an effective field, which is similar to the 
Lorentz local field. The latter is only defined for induced dipole moments, while the former 
is introduced to deal with permanent dipole moments. 

The present consideration can be extended to ferrofluids, which are a suspension [28] 
containing ferromagnetic particles embedded in a carrier liquid. However, the particles in 
a ferrofluid possess a much smaller size than in the ERMR solid. Also, the particles can 
form chains inside ferrofluids under the presence of a moderate magnetic field, without the 
need of a strong magnetic field. For ERMR solids, a strong magnetic field has to be used to 
induce the formation of particle chains. 

Finally, since the particles in ERMR solids are located very close, it is instructive to 
take into account the effect of multipolar interactions between the particles [29-31] on the 
nonlinear ac responses. 

To sum up, we have applied a Langevin model to investigate the nonlinear ac responses 
of ERMR solids. For the longitudinal field case, it has been shown that both even- and odd- 
order harmonics are induced to occur. In contrast, only the odd-order harmonics appear 
for the transverse field case. Moreover, these harmonics can be affected by the degree of 
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anisotropy of the ERMR solid, as well as the field frequency. In particular, the second-order 
harmonics are of several orders of magnitude larger than the corresponding third-order. 
Thus, it is possible to real-time monitor the structural transformation of ERMR solids by 
detecting the nonlinear ac reponses. 
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FIGURES 

FIG. 1. Fundamental, second- and third- order harmonics of Frohlich field as a function of the 
field frequency, for various gL- Parameter: T = 298 K. 

FIG. 2. Fundamental, second- and third- order harmonics of the orientational magnetization 
as a function of the field frequency, for various gL- Parameter: T = 298 K. 

FIG. 3. Same as Fig. 1, but for different temperatures. Parameter: gi = 1/5. 

FIG. 4. Same as Fig. 2, but for different temperatures. Parameter: gL = 1/5. 
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Fig. 2/Huang and Yu 
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